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Abstract 


In turbulence research and flow applications, turbulence models like RaNS 
(Reynolds averaged Navier-Stokes) models and LES (Large Eddy Simulation) 
are used. Both models Alter the governing flow equations. Thus a scale sep¬ 
aration approach is introduced for modeling purposes with the large scales 
simulated using a numerical scheme while smaller scales are assumed to be 
less important and might be modeled more or less easily. Unfortunately 
small scales are frequently of big importance, e.g. in reactive flows, wall 
bounded flows, or flows with significant Prandtl or Schmidt number effects. 
Recent alternatives to these standard models are the class of models based 
on the one-dimensional turbulence (ODT) idea, like ODTLES. The ability 
of ODT to capture highly turbulent flows (recently up to = 6 x 10^) 
allows ODTLES to realize 3D resolutions basically independent of the tur¬ 
bulent intensity. In two papers we provide a formal theory and application 
of an innovative modeling strategy for highly turbulent flows in domains of 
moderate complexity: In part I (see Glawe et ah (2015)) a new general Al¬ 


tering approach, called XLES (extended LES), is introduced. Contrary to 
LES, XLES is based on 2D Altering of the governing equations, whereby ad¬ 
ditional small scale terms are interpreted numerically. In this work a new 
ansatz for the ODTLES model is introduced as one special approach in the 
XLES family of models by incorporating the ODT model into XLES. The 
ODT model introduces microstructures not captured by the XLES Altered 
equations. To illustrate the ODTLES model capabilities, turbulent channel 
and duct flows up to friction Reynolds number Rcr = 10000 are studied. 
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1. Introduction 


In fluid mechanics an increasing number of scientific and industrial prob¬ 
lems are amenable to computer simulations due to growing computational 
power. Realistic problems, e.g. in engineering and meteorology, require the 
application of turbulence models to reduce the computational effort. For 
an overview of turbulence properties and model approaches see e.g. [Pope 


( 2000 ) 


Large Eddy Simulation (LES, see e.g. 


Sagaut (2006)) is widely used in 
fundamental research and increasingly in industrial applications. In LES, 
3D spatially Altered equations, containing large scale properties down to the 
inertial range of the turbulent cascade, are solved numerically. The effect of 
the unresolved scales is modeled, e.g. by an eddy viscosity model (see e.g. 


Germano et ah (1991)). 

In the area of reactive and wall bounded flows, especially with significant 
Prandtl or Schmidt number effects, small scale properties are very important, 
but not adequately represented by simple sub-grid models. 

Recent alternative models including the one-dimensional turbulence (ODT) 
model (see e.g. Kerstein (1999) and Kerstein et al. (2001|) describe the 3D 


turbulence in a ID sub-domain, whereby the numerical representation of 
molecular diffusive effects becomes computationally feasible also in highly 
turbulent flows. E.g. Meiselbach (2015b) presents wall bounded flows up to 


Rcr = 6 X 10^. This model appropriately describes individual flows involv¬ 
ing one characteristic and predominant direction including the full turbulent 
cascade, which is valid inter aha in important problems within the research 
fields of fundamental combustion or atmospheric science. 

In ODTLES the ODT model represents the microstructure terms (of¬ 
ten called sub-grid model SGM), whereby the macrostructural model has to 
fulfill special requirements, which are satisfied by the extended Large Eddy 


Simulation (XLES) ansatz, introduced in part I (Glawe et al. (2015)): Three 


2D Altered and coupled sets of equations are solved and discretized by three 
XLES-grids. Each set of equations corresponds to one predominant Gartesian 
direction. The XLES unresolved sub-grid scale (SGS) terms are connected to 
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the ODT advancement, where ODT is interpreted in terms of Navier-Stokes 
advection, rather than a stochastic (Monte-Carlo like) model. This so called 
ODTLES model (2D hltered XLES eqnations with ODT microstrnctnre mod¬ 
eling) combines the diffnsive and advective small-scale effects compnted by 
ODT with the ability to discretize 3D domains within a macrostrnctnral 
model. Becanse ODT is able to describe the fnll tnrbnlent spectrnm in a 
ID snb-domain, only 3D effects not represented by ODT (e.g. the domain, 
secondary instabilities, ...) need to be resolved in 3D by XLES. 

Here we distingnish between the expressions XLES to describe the XLES 
approach inclnding an approximation or model for arising microscale terms, 
ODTLES if these terms are in particnlar described by ODT and XLES-U 
(XLES nnclosed), if microscale terms are neglected. 

A previons version of ODTLES, introdnced and examined by [Schmidt 
“(|2008|, 


et al, 


Gonzalez-Jnez et al. (2011), and Glawe et al. (2013), solves 


weakly conpled XLES eqnations that generate inconsistent 3D large scale 
velocity helds on the three XLES-grids. This introdnces certain oscillations 
of the root mean sqnare velocity (reported by Gonzalez-Jnez et al. (2011)) 
bnt no fnrther major qnalitative change in resnlts. 

Thongh primarily intended to explain ODTLES to the LES commnnity, 
XLES also gives new insights into ODTLES (e.g. allows to represent scalar 
properties by multiple XLES-grids consistently) and moreover is a novel au¬ 
tonomous modeling strategy because the ansatz is very general and not lim¬ 
ited to one-dimensional models like ODT. 

There are other approaches like LES-ODT (e.g. by Cao and Echekki 


(2008)), LES-LEM (by Menon and Kerstein (2011)), and LEM3D (e.g. by 


Sannan et al. (2013)) connecting ID turbulence models and 3D computations. 


These approaches are not considered in detail in this work. 

In this work the ODT model is briefly introduced in section ^ and the 
XLES approach, introduced in detail in part I (Glawe et ah (2015)), is sum¬ 
marized in section [H Section |4] introduces the ODT model as closure for 
XLES. The time advancement cycle of the ODT closed XLES model (OD¬ 
TLES) in section [^summarizes the model derivation, followed by channel and 
duct flow results in section [^ comparative estimations of the computational 
costs in section and hnal conclusions in section [^ 

In accord with part I ( [Glawe et al. (2015)) we assume the flow to be de¬ 
scribed by the incompressible Navier-Stokes equations for a Newtonian fluid 
with constant kinematic viscosity {u) and constant density (for simplicity we 
assume po = 1). 
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Note that in this work no Einstein summation convention is used. 


2. One-Dimensional Turbulence Model (ODT) 


The ODT model describes the dynamics of a three dimensional turbulent 
flow within a one-dimensional sub-domain, including fully resolved molecular 
diffusion. Thus ODT is a dynamical model, able to describe e.g. the turbu¬ 
lent channel flow including high order flow statistics with high accuracy (see 


Appendix A). 


ODT stand-alone is able to compute meaningful results, even with one ve¬ 
locity component (see Kerstein (1999)). Nevertheless to capture anisotropic 


flow behavior and especially as a closure within a 3D approach, two or three 
velocity components are advantageous. 

In wall-bounded flows (especially in turbulent channel flows) described 
by ODT with three velocity components, the wall-normal and the spanwise 
velocities are identical, so two-component ODT captures similar statistical 
flow properties. (Note that ODT results in Appendix A are computed with 
2 velocity components). 

In this section the ODT time advancement is described briefly. Here we 
introduce ODT including 2 velocity components. This ODT model is a mod- 


ihcation of the ODT vector formulation by Kerstein et ah (2001) including 
3 velocity components. 

ODT emulates the time evolution of a turbulent 3D fluid in a ID sub¬ 
space, which is oriented in the Cartesian x^-direction. The 2 ODT velocity 
components Uk,i (with k ^ i) are oriented orthogonally to the x^-direction. 
The time evolution of a velocity held Uk,i in this ID subspace is described 
by: 


{dtUk^i 'DoDTk(wfc,i) + ek,i{uk,h Xq, 0) = 0 with k,i = {1,2,3} A i ^ k (1) 

with the ODT diffusion term 'DoDTk('^fc,i) = ~^dl^Uk^i which is numerically 
approximated by an implicit Euler scheme in time and a central difference 
scheme in space. The index notation resembles the XLES index notation used 
in the following sections with the velocity Uk,i oriented in the Xj-direction 
within a ID sub-space oriented in the x^-direction {i ^ k). 

The term ek^i{uk/,XQ,l) is an instantaneous eddy function affecting Uk^i 
within the eddy range Xk G [xo,Xo + 1]- The maximum eddy length is 
enforced, hence I < The eddy function Ck^i is introduced to represent a 
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stochastic procedure that emulates turbulent advection: 

ek,i ■ Uk,i{xk,t) -)■ Uk,i{f{xk,l),t) + CiK{xk). ( 2 ) 

Note that the advection function Ck^i depends on both velocity components 
Uk,i (with i = {1,2,3} A z 7^ k) due to q in Eq. (|^. The mapping func¬ 
tion f{xk,l), representing fluid transport, is measure preserving (the non¬ 
local analog of vanishing velocity divergence), continuous, and satishes the 
requirement of scale locality (at most order-unity changes in property gra¬ 
dients). These indispensable physical requirements for f{xk,l) are satished 
by a triplet map, which places three compressed copies of the original prohle 
{uk,i{xk), Xk € [xq^Xq + /]} in the eddy range. The middle copy is reversed 
to preserve continuity. The triplet map /(x^, 1) ^ Xk is: 


f{Xk,l) 


-■ xo + 


3{xk - Xo) 

{Xk - Xo), 


( 3 ) 


else. 


In Eq. (|^, K(xk) is a kernel function which in combination with the 
amplitudes Cj assures momentum and energy conservation and controls the 
energy redistribution among the velocity components. A possible dehnition 
is: K{xk) = Xk — f{xk,l). This energy redistribution is a ID interpretation 
of the pressure-fluctuation effect in a 3D flow and therefore is called pressure 
scrambling. 

Determination of the amplitudes Cj requires additional modeling. 


Kerstein et ah (2001) derive the amplitudes: 


27 


D = ^ I -UK-,k,i + sign(Mi^;fc,i) | ^ ^ k, j ^ k 


with the dehnition 


UK-,k,i =p Uk,i{f{xk))K{xk)dxk 


( 4 ) 

( 5 ) 


and the transfer matrix 


aT = a 


-1 1 
1 -ir 
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( 6 ) 








The free parameter a ensures the amplitudes c* in Eq. Q to be real values 
for 0 < a < 1. We choose a = 1/2 corresponding to the equalization of the 
two component available energies in the present formulation. 

During an eddy event, the transfer matrix T redistributes the turbulent 
kinetic energy among velocity components. 

This ‘pressure scrambling’ accounts for the tendency for pressure fluc¬ 
tuations to restore isotropy and is invariant under exchange of indices. By 
construction the momentum and total energy are not changed by the pressure 
scrambling. 

During the ODT time evolution in Eq. Q, the eddy size I and the location 
Xq are sampled from a probability distribution representing the physics. For 
a given {l,xo} an eddy turnover time Xg can be calculated leading to an 
occurrence frequency —. Since the ODT triplet map is an instantaneous 
process, the frequency for the eddy specihed by {l,xo} is chosen from an 
event rate distribution: 


C 


C 


PTe{Xo,l) 


(7) 


involving particular definitions of the turbulent kinetic energy E^in and the 
potential energy Epot- The latter vanishes for the cases considered in this 
work. Note that the work by Ashurst and Kerstein (2005) introduces a 
variable density formulation of ODT including Epot- The values C, 

and Z are adjustable model parameters. The latter is introduced to cut off 
eddies with unphysically small energy and the parameter C is an overall rate 
coefficient determining the strength of the turbulence. The maximum eddy 
length is chosen to characterize the largest (global) scale within the flow, 
e.g. the channel half height. 


3. Extended Large Eddy Simulation (XLES) 


In sections 3.1, 3.2, and 3.4 the XLES framework is summarized briefly. 


to allow focus on features related to ODTLES without repeating details in¬ 
troduced in part I (Glawe et ah (2015)). 


Additionally a time scale separation is introduced in section |3.3| to sim¬ 
plify the interpretation of microscale terms with reference to the ODT ad¬ 
vancement. 
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(a) [hhhjui = (b) grid 1: for ui^i. (c) grid 2: for M 2 ,j. (d) grid 3: for 


Figure 1: In XLES the velocity components Ui are resolved using multiple 
XLES-grids illustrated in [Tb||Td| 3D large scale properties, corresponding to 


a standard LES grid with (illustrative) X'les = 4 cells, are illustrated in la 
Illustrative XLES-resolved small scale (‘RSS’) properties are approximated 


using A^rss = 16 cells in lb Id 


3.1. XLES: Spatial Filtering 

The basic XLES concept is to apply 2D Liters to the velocity held, main¬ 
taining one Cartesian direction highly resolved. Using these filters, three 
2D hltered velocity helds, each corresponding to one highly resolved Carte¬ 
sian direction, are introduced and discretely represented by three staggered 
XLES-grids, illustrated in hgure [Tb] - [Td} 


A vector-matrix notation (which we will refer to as XLES-vector notation) 
is used to represent the three 2D hltered velocity helds: 


j 

[hl 2 ]UiJ 


fhh 0 ^ \ ( “A 

j 0 I 1 I 3 0 j Ei = Ui, with, i = {1,2,3} (8) 

V 0 0 hhj \uij 


where i indexes the velocity components. We introduce another index k 
referring to the velocity held Uk,i represented in XLES-grid k (similar to the 
index notation in section]^. 

Here without loss of generality (w.l.o.g.) the 2D hlter operator [Z 2 / 3 ] (cor¬ 
responding to XLES-grid 1) represents a tensor product of ID hlter operators 
[ 12 ] and [I 3 ]. This tensor product is commutable: [/ 2 ^ 3 ] = [^ 3 ^ 2 ]- 

The ID hlter operators are used to derive micro-and macroscale velocity 
terms: 


Ui — [hhh + S 1 I 2 I 3 + ^ 1 ^ 2/3 + ^ 1 ^ 253 ] Ui (9) 

+ [-5152/3 + 51/2S3 -|- /1S2S3 + S1S2S3] Ui . 

"-V-' 

=S (denoting XLES SGS terms) 
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with the ID small scale operator = 1 — (A; = {1,2,3}) and the unity 

operator 1. 

These scales are decomposed (using the 2D hlter matrix 


( [hh]ui\ / [hs2h + hhs3\ui\ 

[hh]'^i j + I [51^2^3 + hhS3]'<J‘i j 

[hh]Ui/ \[si/ 2^3 + 
into three terms: 


/[5]mA 

+ = f ^ u, + + [5] 

\[S]uil 

( 10 ) 


1. ‘Directly Resolved’: 

Mj is discretely represented in the XLES-grids. 

2. ‘Indirectly Resolved’: 

The resolved small scale (RSS) velocities 

h, ^ = h, - ( 11 ) 


(with = [/i/ 2 ^ 3 ]Mj) lead to coupling terms between the XLES- 

grids, and are directly resolved by another XLES-grid (see section 3.2 


for details). The coupling matrix C and the small scale matrix are: 


i 

(0 

1 

1^ 

\ i 

(Sl 

0 


C = 

1 

0 

1 

and = 

0 

S2 

0 

— 1 

u 

1 

0) 

1 “ 1 

^0 

0 



( 12 ) 


The LES velocity held corresponds to the ID hltered XLES ve¬ 
locity held: [4]hfc,i (see hgure [^. 

3. ‘Not Resolved’: 


[S]Ui = [S 1 S 2/3 + S 1 I 2 S 3 / 1 S 2 S 3 -f SiS 2 S 3 ]Mi = Ui (13) 

is not resolved in any XLES-grid and leads to the XLES microscale 
terms. 

Because the microscale model (ODT) is able to represent large scale ef¬ 
fects and the macroscale model (XLES) contains ID small scale terms, the 
classical term ‘scale separation’ is misleading for XLES and especially OD- 
TLES. Thus we will refer to a ‘hlter separation’. 



3.2. XLES: Momentum Conservation 

The 2D filtered XLES momentum equations are 


0 + I d, 



.u. * u. 

3 —J —I 


( 14 ) 




XLES 


i=i 

I p^XLES^ _j_ ^spatial 


i=i 

with the entry-wise multiplication * between XLES-grid vectors (and ma¬ 
trices). The decomposed SGS Reynolds stresses and the cross-stress 

terms are combined to obtain the XLES sub-grid scale (SGS) advection 

terms: 

f + + UiUj) ^ (15) 

and its coupling * C Here Uj is the vector notation for Eq. (13), 

hj is dehned in Eq. (11), and l) is dehned in Eq. (17). 

ing stress terms of the resolved small scale velocities 

(16) 


The linearized coup 


which has a discrete representation. Hereby the matrix 1} with 


= cHi 


K% ii% 


1 






2 '■S 
3 ^^2 1 


(17) 


includes the deconvolution operator which reconstructs hltered informa¬ 
tion. Due to the coupling the large scale influenced by highly resolved effects 
in a specihc XLES-grid is communicated to the other XLES-grids. Hereby 
the reconstruction is executable because the operator is only applied to 
large scale terms in x^-direction. For details about the coupling procedure 


and an possible algorithm for a discrete deconvolution, please see (Glawe 


et al., 2015, section 3.1). 


The spatial XLES model error terms are: 

^ dxjl_ jl(hlj'U2,i T T T T ^3,ji^l,i T ^3,j^2,i) T C^ij 

( 18 ) 


including non-linear coupling terms and the 2D Leonard stresses 


2D 
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3.3. XLES: Time Scale Separation 

The XLES advection terms are represented by three overlapping XLES- 
grids, including coupling terms between these XLES-grids, and an additional 
ODT advancement for ODTLES. This ODT advancement involves instanta¬ 
neous stochastic mappings whose instantaneous nature is in conflict with the 
idea of high order time integration schemes. 

On the one hand, the simplest and physically most convenient way to ad¬ 
vance the coupling terms and the dynamical SGM within an XLES framework 
is an explicit Euler scheme, which allows a straightforward interpretation of 
coupling terms and the stochastic turbulent advection within ODT. 

On the other hand, an efficient numerical advection scheme includes high 
order time integration, which is even required for stability reasons by some 
spatial discretizations. 

A known compromise is to linearize the advection; The linear advection 


part is advanced by a high order numerical scheme (details in part I (Glawe 


et al. (2015))), while the non-linear part is implemented by a 1st order explicit 


Euler scheme. One possible way to interpret such an approach is to integrate 
the dynamical velocity held over one time step and use this velocity held to 
advect the dynamical variables within the next time step. This approach 
can easily include random occurrences of instantaneous ODT mappings and 
therefore is suitable for ODTLES. 


In part I ((Glawe et ah, 2015, section 3.3)), alternative time schemes are 


suggested which potentially avoid the linearization of the advection terms. 

In contrast to RaNS models, time averaging is not applied to the dy¬ 
namical variables, but to the advecting velocity variables. This is especially 
reasonable because ODT is a dynamical model introducing small time scale 
ehects. 

Formally a time scale separation is invoked that simplihes the inclusion 
of ODT into XLES, though this procedure is not required by XLES (in part 
I (Glawe et al. (2015))) time scales are not separated). The separated time 
scales are: 


Uj = {Uj) + {Uj} (19) 

with the large time scale 

1 

(uj) = - Uj dT ( 20 ) 

^ Jt 
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and the small time scale (fluctuations) {uj} = Uj — (uj). Note that the time 
Alter and spatial Alters in XLES are independent of each other:([Zfc]Mj) = 
[Ik] (Uj)- 

The integral time r corresponds to the 3D large scale flow. A natural 
choice for r is the discrete time step size of the 3D large scale advancement 
scheme. Thus the modeling strategy is directly connected to the numerical 
realization. 

If an explicit time integration scheme is used, r is numerically restricted 
by a (global) Courant-Friedrichs-Lewy condition. Here CFL defines the con¬ 
stant CFL-number depending on the implemented numerical scheme: 

r = CFL min ((21) 

within all XLES-grids k = {1,2,3} and with all velocity directions Xi {i = 
{1,2,3}. Local time step restrictions are possible but not implemented. 

The time scale separation within XLES-U implies additional error terms 
caused by time scale separation: 

= 1:4, ( 22 ) 


By defining the integral time scale r based on the resolved small scale cell 
size in XLES-grid k, the time scale separation is suppressed, because 

all XLES-U velocities are large scale in time: 


r = CFL min 

k,i 




(23) 


The UFL-number can be used to switch between Eq. (^21h and Eq. (23) and 


thus becomes a model parameter balancing (and controlling) the temporal 
model error model performance. 

Thus the UFL-number can be increased for performance reasons: OD- 
TLES is still stable and well defined for CFL = 0.5^)^ which e.g. for 
A^rss = 512 and Ales = 16 leads to a factor 32 increased XLES time-step 
size (ODT advancement is indirectly influenced) with additional model errors 
(see Eq. (22)) having only a small impact on specific problems. In this work 


all results are computed with considerably small UFL-number, following Eq. 


(23). 
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Since ODTLES allows a huge number of turbulent ODT events within 
the time scale r, the averaged velocities are smoothed, which is part of the 
ODTLES modeling strategy. Again the OFL-number is controlling this mod¬ 


eling impact, because with r based on Eq. (23) the number of ODT turbulent 
events is decreased. 

ODT describes fluctuations (small time scale terms) in ODT-direction Xk 
corresponding to advection terms of the form dj{uf, k}^k,i- The corresponding 
time averaged advecting velocity {uk,k) is specihed due to mass conservation, 
see section 3.4 The terms correspond formally to the spatial 

XLES macroscale (they are part of the model error in Eq. (22), controlled 
by CFL), but can be interpreted by ODT, which is possible due to the time 
scale separation. 

Additionally to the terms ODT directly and indirectly (due to cou¬ 
pling) represents fluctuations across 3D large scale cells: 




9xi {^ 1,1 }^ 1 ,! 

{dxjiVLj} * + {QpXjiikj} * = {g + ij * Q) I dx2{u2,2}'^2,- 

dx3{u3,3}u3,-, 


I temporal 
+ ^ODTLES- 


(24) 


This model assumption leads to a model error term replacing Eq. (22): 


^ dx2{u\^2]u\^i + dx^{ui^3}ui^i 

^ODTLES ~ 1 + dx^{u2^3\u2^i 

^dxA'^3,l}u3,i + dx2{u3,2}u3,i 


(25) 


The error term in Eq. (25) summarizes all fluctuating terms not in ODT 
direction Xk- Again this error term is controlled by the OFL-number. 

Applying the time scale separation to the advecting velocities leads to a 
modihed XLES momentum equation (compare to Eq. 0): 


0 =d. 


+M 





* u. 


J =1 


ODT + A 


a 


temporal 


(26) 


((«_,)» a-(«,.) 

i=i 




+ lUCM 


^iXEODT- 
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The advection terms assumed to be modeled directly by ODT (involving 3 
velocity components) are: 


M 


ODT 


+ d^AP} + * pLi) 

i=i 


(27) 


The derivation of JSAodt is tailored for ODT, which is emphasized by the 
acronym ODT. Additionally the ODT advancement is coupled between the 
XLES-grids ([i *£A1 odt)- 

Note that spatial large scale terms of the form included 

in the ODT coupling P * C 


The term dxAP} = L^^xAp} describes pressure fluctuations. By con¬ 
struction ODT is mass conservative, nevertheless pressure fluctuations are 
modeled by applying the so called pressure scrambling (see section]^. 

The full ODTLES advancement cycle is summarized in section 
An alternative approach, applied by Cline (2015) within the lattice-based 
multiscale simulation model (LBMS), is to couple each individual turbulent 
event within the ODT advancement, instead of the time averaging approach 
introduced here. This approach potentially introduces small time scale com¬ 
munication within a parallel algorithm. 


3.4- XLES: Mass Conservation 

In the incompressible flow regime, the 2D hltered velocity helds need to 
be divergence free to ensure mass conservation. 

Because the XLES dynamics take place on the integral time scale (r), the 
2D hltered mass equation 

0 = X = X S-. (5)“" + X (a) (28) 

2=1 2=1 2=1 


is enforced for the integral time scale by the procedure described in this 
section, while velocity huctuations (corresponding to the small time scale), 
described by the ODT advancement, are mass conservative by construction 
(see section]^. 

The mass conservation of the 3D large scale velocity is enforced 

by a standard approach; A pressure Poisson equation is solved, leading to 

a large scale pressure held (p) (see standard textbooks, e.g. Ferziger 


T3 








and Peric (1999)). The resulting pressure gradient {p 
divergence free velocity field by solving: dt 


LES 


enforces a 

/ \ LES 

+ dxi = 0 

(see section]^ and especially Eq. ([^). Here a consistent 3D large scale field 
= [/i] {ui,i) = [h] {u 2 ,i) = [^ 3 ] {u 3 ,i) is required, which is enforced by 
the coupling terms and * C AIqpt)- 

The resolved small scale (RSS) divergence X]i=i (Uj) = 0 vanishes if a 


2 

2 


{ui) dx'j^ is used (proof in (( 


Glawe et ah 


2015 


box hlter [h] (ui) = ^ J 
Appendix B))). 

Since the 3D velocities are divergence free (after solving the pres¬ 

sure Poisson problem), a direct solver can compute one velocity component 
{uk,k) in each XLES-grid k within one 3D cell of the size (w.l.o.g. in 
XLES-grid 1): 




-Axi 


+ Xi = (mi) 


- \LES 


/S.Xi 


(29) 


-+Xl 


Ax-i 

■ ~T~ 


dx2 (“1,2) dx) - 


-+Xl 


■ ~T~ 


dx3 (hi, 3 ) dx[ 


for Xi < Axi- Owing to the absence of a 3D small scale velocity held, small 
scale pressure effects vanish from the equations: (p) = but ODT 

explicitly models small time scale pressure effects (see section]^. 

Because 9 velocity components within 3 XLES-grids are redundant, one 
dynamical velocity component in each XLES-grid is omitted by multiplying 
a matrix of Kronecker deltas 

A -5u 0 0 

0 1 - ^ 2 * 0 
0 0 1 - 


1 - = 


with (1 - 5ki) = 



(30) 


to the XLES momentum equation Eq. (26), leading to 

0=(l-ij)4,(p) +{di-v'^dl\(t-y)u. 


(31) 


i=i 


+ (| - 4) 11 {^-1%) * f * dx,g {{%) * Ui - 




-p A^ Qprp -\- * C -|- (l - ^.)P- 
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Figure 2: Illustrative ID stack in XLES-grid 2 in a; 2 -direction at {xf,x§) 


Here terms to be modeled by the ODT model are 

Modt = (1 - i.) E + 4.{p} + d,.{u,} » «,) (32) 

and involve two velocity components orthogonal to the highly resolved direc¬ 
tion in each XLES-grid (introduced in section]^. 


4. From XLES to ODTLES 


In this section ODT is interpreted as a microscale model within the XLES 
approach. In ODT a stochastic process mimics 3D turbulent advection within 
a ID sub-domain. Thus, an interpretation of ODT in terms of the Navier- 
Stokes advection is not straightforward and ODTLES is not directly de- 
ducible from Navier-Stokes equations. 

Nevertheless XLES microscale terms w.l.o.g. in XLES-grid 2 odt 


Eq. (27)) are interpreted by IVlesi ^ -^less so called stacks, each containing 


highly resolved ID information, e.g. defined by one line at {xf,X 2 ,x§) with 
constant xf and xf (see hgurej^. The microscale terms A^ 2 ,onT(^i) ^ 2 , ) 

in each of these stacks S contain 2 velocity components and can be modeled 
by the ODT advancement (Eq. 


^2, = e2,i{u2,i] X2,1) + (JoDTk ) ^nd I < (33) 


with the ODT model error (JoDTk- Note that the ODT model introduced in 
section also advances 2 velocity components which are orthogonal to the 
Xfc-direction. 
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The maximum eddy length within the ODT model corresponds to 
the largest (global) scale. By using ODT as a model within the XLES frame¬ 
work, the maximum eddy size dehnes the boundary between turbulent 
scales described by the 3D advection scheme and by the ODT turbulent ad- 
vection. Since ODT should capture turbulent effects not resolved by the 3D 
advection scheme, the maximum eddy size mainly depends on the nu¬ 
merical properties of the 3D advection scheme and needs to be determined by 
numerical tests (not shown here). We found to be convenient, 

which corresponds to the ability of the implemented numerical 3D advection 
scheme (see (Glawe et al., 2015, section 3.3)) to resolve e.g. the Kolmogorov 
length scale with approximately one 3D cell in the ‘DNS-limit’ of XLES (all 
scales are represented in 3D) as shown in part I ((Glawe et ah, 2015, section 
3.2)). 

Glosure of XLES can involve any form of modeling that specihes the RSS 
time advancement on an entire XLES-grid such as XLES-grid 2 shown in 
hgure 2. This is not required to involve a collection of model instantiations on 
individual stacks, such as the illustrative stack in that hgure. Nevertheless, 
sub-grid ODT within ODTLES is formulated in this way. On this basis, 
the XLES 3D advection can be viewed as a form of coupling of the ODT 
instantiations within one grid. In this context, the grid-to-grid coupling can 
be seen as a higher level of coupling. This is mentioned because previous 
ODTLES formulations did not envision the XLES framework with ODT not 
being the only conceivable RSS closure strategy within an XLES grid 

In previous ODTLES formulations, pressure projection was the only grid- 
to-grid coupling. The insufficiency of this coupling is illustrated by a notional 
extension of the method to include passive scalar properties. These are not 
subject to pressure projection, so passive scalars on different grids would 
not be coupled at all, an obviously unsatisfactory situation. The present 
approach is straightforwardly extended to scalar properties in a manner that 
provides appropriate grid-to-grid coupling. This conceptual advantage of 
XLES does not imply that previous ODTLES formulations, which do not 
include scalars, are necessarily dehcient in some way. Indeed, for reasons 


that are not entirely clear at present, they appear to perform well (Schmidt 


et al. (2008), Gonzalez-Juez et ah (2011)), with the caveat that the latter 


documented oscillations in velocity root mean squares that are not seen in 
results obtained using the present formulation (see sections 6.1 and 6.2). 
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4.1. ODTLES: ODT-limit 

ODTLES includes another distinguished limit that we refer to as the 
‘ODT-limit’: On the one hand ODTLES collapses to the ODT stand-alone 
model if only one 3D large scale cell represents the full domain. On the 
other hand the XLES microscale terms correspond to the full Navier-Stokes 
equations in this limit. 

Thus the ODT model error (Todt can be estimated by comparing ODT 
and DNS results for a turbulent channel (see Appendix A). 

The ability of ODT to describe the full spectrum of 3D turbulent effects 
is a required property to get an ODTLES 3D resolution largely independent 
of the turbulent intensity unless Reynolds-number variations trigger a global 
flow structure transition (see duct flow in section 6.3). Indeed, demonstrated 
model performance in the ‘ODT-limit’ (i.e. ODT stand-alone, see Appendix 


0 strongly indicates that ODT adequately describes the XLES model terms 
This is also supported by ODTLES results that are shown to be in good 
agreement with DNS in sections |6.1[ |6.2[ and |6.3[ Unfortunately a detailed 


theoretical investigation of the ‘ODT-limit’ requires a convenient ODT in¬ 
terpretation in Navier-Stokes terms, which is not derived to a satisfying level 
yet, but ensemble statistics are formally analogous to corresponding Navier- 
Stokes terms to a considerable extent (e.g. the interpretation of ODT budget 
terms of the turbulent kinetic energy in Appendix A). 


4.2. ODTLES: ODT Modeling Effects 

ODT introduces local turbulent events depending on the local flow state. 
In low Reynolds number channel flows the 3D grid is under-resolved only in 
the near-wall region (unless the grid is very coarse) and thus ODT works as a 
dynamical and highly accurate near-wall model, as hgure illustrates. For 
highly turbulent flows, the 3D resolution in the core region of the channel 
is under-resolved too: In this case ODT small scale eddy events additionally 
occur in the core region introducing local turbulent transport effects, as hgure 
l3bl illustrates. 

In contrast to wall-modeled LES, the ODT turbulent transport treats 
all regions consistently without introducing additional assumptions for the 
near-wall region. 

Additionally the ODT modeling depends on a fully resolved (ID) how 
state and thus allows dominant small scale ehects (e.g. local stratihcation, 
chemical processes, ...), which are not well captured by commonly applied 
eddy viscosity models. 
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V(v/u,^) 

(a) Eddies at Re-r = 395. 


t/(v/u,") 

(b) Eddies at Rcr = 2040. 


Figure 3: ODT turbulent events (‘eddies’) occurring in XLES-grid 2 (super¬ 
position of eddies in all ODT domains). The size of the 3D cell closest to 
the walls is illustrated by red lines (Xles = 32). For low Reynolds numbers, 
ODT primarily acts as a near-wall model (3a); for high Reynolds numbers 
ODT acts as a sub-grid model over a larger extent of the flow domain (3b). 


4.3. ODTLES: Coupled ODT Models 

The turbulent ODT advection Ck^i (in XLES-grid k) models the unresolved 
XLES terms These are coupled across the XLES-grids due to the 


SGS coupling term p * O AIqdt (^^e Eq. (31)). The ODT diffusion terms 
’^ODTk = (for k i) m. XLER^id k represent the molecular 

diffusion as a continuum and are connected to the XLES diffusion terms 
(1 — 5^)(—zz) w.l.o.g. the XLES diffusion terms in grid 1 


ydl ui^i + + I'dlJii.i , withi ^ {2,3} 


(34) 


are solved in three ways (compare to section 3.1): 


1. Diffusion directly resolved by ODT: 

The terms with i = {2,3} are interpreted by ODT incorpo¬ 

rated in XLES-grid 1. These terms are resolved by Xrss cells (repre¬ 
senting molecular diffusion, similar to DNS) 

2. Diffusion indirectly resolved by ODT: P * ODqdt 

The terms —are interpreted by ODT domains 
residing in XLES-grid 2 respectively XLES-grid 3. W.l.o.g. the first 
term is coupled from XLES-grid 2 to XLES-grid 1 by — 

(index notation for P * O^odt)- The diffusion is fully resolved, but 
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additionally filtered (convolved) and deconvolved (see (Glawe et al. 


2015, section 3.1)). 


3. Diffusion resolved by XLES: ^xles 

The terms —vdl^Ui ^2 and are not interpreted by ODT in 

any XLES-grid (this is caused by (1 — 5.); see section 3.4). A nu¬ 
merical interpretation is possible within XLES-grid 1 using IVles cells: 

with i = {2, 3}. These diffusive terms are not resolved down 
to the molecular level. The XLES resolved diffusion terms are written 
as ^xLES numerically represented by an explicit Euler scheme in 
time and a spatial central difference scheme. 

In summary the ODT model is incorporated into XLES-grid k by inter¬ 
preting diffusive effects "DooTk and the microscale advection terms AdoDT^- 


M 


ODTk 


- i^dl Uk,i ^ ek,i + VoDT{uk,i)y fori 7 ^ k 


(35) 


Additionally the (diffusive and advective) ODT terms are coupled between 
the XLES-grids by * Q{ei{iki) + 2oDT(ili))- XLES diffusion terms Vxles 
are introduced to represent diffusive terms not captured by ODT. 

The under-resolved diffusion terms ^xles generally much smaller than 
the correct local diffusion and might be omitted in typical applications. Nev¬ 
ertheless these terms are conceptually desirable, because they allow the cor¬ 
rect behavior in the ‘DNS-limit’. 


5. ODTLES: Time Advancement 

To advance the XLES equations in time a modified predictor-corrector 
procedure is used: The XLES momentum equations including coupling terms 
and ODT advancement are solved, predicting velocity helds in each XLES- 
grid. Simultaneously time averaged velocity helds are computed. A corrector 
step enforces the time averaged velocity helds to be divergence free (to ensure 
mass conservation). 

Since the predictor-step involves ODT advancement and several coupling 
terms, a fractional time step algorithm is introduced: 
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1. predictor: 


Ui=Ui{t) + j^ d^^(p^^dt' 

/*t+T 3 

dt' 

■'t j=i 

rt+T 3 

+ / [U ^ a,, - (m ^ 

,=1 ^ 


** /v >(: I 

“i =“i + 


a,’” =ar 


(36) 

(37) 

!7 !7.*.I-ES'\ 




(38) 


r^t+T 


r^t+T 


fi."” =*(” + / (eitef) + ®odt) (fif) d«' + / VxLEslm") *' 

(39) 

rt+T 

Ut + r) =ir" + / l^*C {ei{u***) + Vodt) {u***) dt' 


2 . corrector: 

ID filter 

solve 


I— \ijlis jlJJ /« \ 

\lkUt+T \—i)t+T 

3 

0 = 5^, (“*)t+r ^ (p 


0 0 


with (p) = 0 I 


0 0 /o 


Z =1 


LES 


LES 


t+T 


(40) 


(41) 

(42) 


2 ^ 0 (p) (no small scale pressure field 

n 7-1 / 


within XLES). 


The subscript (e.g. () J introduced in time averaged properties indicates the 
time t (averaged over the last time step) respective t + t (averaged over the 
actual time step with the time step size r). 

In the ODTLES advancement cycle the ODT advancement is the most 
costly sub-process which leads to a highly parallelizable algorithm. 


The concrete numerical implementation of the individual steps Eq. (36 


Eq. (40) is shown in Appendix B 
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6. Application: Channel Case 


6.1. ODTLES: Convergence Study 

To verify the ODTLES model, we conduct a numerical convergence study 
by computing a turbulent channel flow with friction Reynolds number Rcr = 
395. 

In part I ( Glawe et al.| 2015, section 3.4) XLES-U (unclosed XLES) 
and LES-U (under-resolved, unclosed LES) results are compared to DNS by 
Kawamura et ah] ( |1999 ) (online available: Kawamura (2014)): both LES-U 
and XLES-U converge towards DNS with increasing 3D resolution, but only 
XLES-U is able to represent the laminar sublayer independent of the 3D 
resolution. (The numerical schemes for LES-U and XLES-U are identical). 

The ODTLES results are compared to XLES-U to show the significant 
effect of ODT as a SGM (see hgure |^. Additionally the DNS results are 
presented for comparison reasons. 

For this convergence study, both XLES and ODTLES are using Ales = 
{16,32,64} equidistant 3D large scale cells and Arss = 512 cells resolving 
the small scale. 

The XLES time step size is limited by Eq. (23) with the Gourant- 


Friedrichs-Lewy number CFL = 0.45 and the small scale cell size 
(the ODT advancement is only indirectly influenced by this GEL condition). 
This choice minimizes the temporal XLES error terms introduced 



length in XLES-grid k equals the 3D large scale cell size 

To produce statistically signihcant results, the flow is averaged for tave > 
25 non-dimensional time units (compared to tave = 20 for DNS) after reaching 
a steady state. 

The mean velocity profiles computed by ODTLES (see figure 4b) and 
XLES (see figure [ia]) are compared to DNS. 

Additionally the streamwise and spanwise velocity RMS (see figure 
and the budget terms of the turbulent kinetic energy (see figure [4d| - [4f[ ) 
shown. 


are 


As derived in Kerstein et al. (2001) ODT pressure fluctuations cannot 
be distinguished from other turbulent transport terms. Thus we combine 
pressure terms into the turbulent transport also in ODTLES. 


21 

















(a) XLES : law of the wall (112,i)- Pro- (b) ODTLES: law of the wall (112,i)- 
files shifted with increasing Profiles shifted with increasing 




(c) Streamwise (ttRMs) arid spanwise (d) Production (Prod) and Dissipation 
(r^RMs) velocity RMS. (Diss) of the turbulent kinetic energy. 




y+ y+ 

(e) Viscous transport of the turbulent ki- (f) Advective transport of the turbulent 
netic energy (tv). kinetic energy (ta). 

Figure 4: Turbulent channel flow results for DNS (small crosses), ODTLES, 
and XLES-U with A^les = 16 (dash-dotted), A^les = 32 (dashed), A^les = 64 
(solid). The small scales are resolved using AIrss = 512 cells. The flow 
statistics are based on the velocity held U 2 ,i- 
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XLES-U shows convergence towards the DNS results with increasing 3D 
resolution. ODTLES is able to represent the flow fleld and turbulent statis¬ 
tics even with the very low 3D resolution iVLES = 16, including the laminar 
sublayer near the walls and the budget terms of the turbulent kinetic energy. 
This significant improvement for channel flow results with very coarse 3D 
resolution compared to XLES-U suggests that ODTLES can encompass a 
wide Reynolds number range for fixed (coarse) 3D resolution. With increas¬ 
ing 3D resolution ODT-speciflc issues in the near-wall statistics decrease in 
ODTLES. 


6.2. ODTLES: High Reynolds Number Flow 

ODTLES combines the ability of the ODT model to describe all scales of 
highly turbulent flows within a ID sub-domain of the full 3D domain with a 
coarse grained XLES approach representing the domain of e.g. a turbulent 
channel and introducing additional 3D effects compared to ODT stand-alone. 

As shown in a convergence study in section 6.1 low Reynolds numbers 
are well described by ODTLES with only IVles = 16 3D cells. To demon¬ 
strate the ODTLES ability to describe highly turbulent flows within a simple 
domain, we conduct turbulent channel flow computations with IVles = 32 
cells (in 3D) and up to N-^ss = 16384 cells to represent additional small scale 
effects, which allows Reynolds numbers Rcr < 10000. 

The CEL number is chosen following Eq. (23) with CFL < 1 and ODT 
parameters are: C = 6.5, Z = 330, and l'^^^ = 

The mean velocity profiles for several friction Reynolds numbers are illus¬ 
trated in figure [5a] and compared to DNS by||Kawamura et ah (1999) and Lee 


and Moser ( |2014a ) (online available: Kawamura (2014) and Lee and Moser 


(2014b)). The ODTLES computation with Rcr = 10000 is in good agree¬ 
ment with the laminar solution near the wall and the law of the wall with a 
von Karman constant k = 0.384, as obtained by Lee and Moser (2014a) for 
Rer = 5200. 

The streamwise and spanwise velocity RMS (see figure 5b) are in good 
agreement with the available DNS in the laminar region near the wall and 
beyond the first 3D cell. The size of the first 3D cell is illustrated by vertical 
lines for the different Rer values. Within the first 3D cell, ODT typically 
has some issues in representing the velocity RMS (this also applies for the 
budget terms of the kinetic energy). 

In figure [5c1 - [5e| the budget terms of the turbulent kinetic energy are shown 
to be in good agreement with the DNS results, especially for highly turbulent 
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(a) Law of the wall. Profiles shifted with (b) Streamwise (urms) and spanwise 
increasing Rct- (—'u^rms) velocity RMS. 




(c) Production (Prod) and Dissipation (d) Viscous transport of the turbulent 


(Diss) of the turbulent kinetic energy. 



kinetic energy (tv). 


. . .DNS, Re 

=1240 

_ ODT, Re^ 

=1020 

— Nles-'S 

Re^=1020 

— Nles-32 

Re^=1020 

— Nles-“ 

Re^=1020 

,-,-Nles-32 

Re^=395 
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(e) Advective transport of the turbulent (f) Wall shear stress tw'- 

kinetic energy (ta). linear (left) and logarithmic (right). 


Figure 5: Turbulent channel flow results for ODTLES (solid) and available 
DNS (dashed) for A^les = 32 (if not indicated otherwise) and with A^rss = 
512 (for Rer = 395), A'rss = 2048 (for Rcr = 1020), IVrss = 4096 (for 
Rcr = 2040), A^rss = 8192 (for Re^ = 5200), and A'rss = 16483 (for Re^ = 
10000). Figure [^compares the ODT^d ODTLES wall shear stress {Rcr = 
{395,1020}) for various IVles to DNS (Rcr = 1240). 












































































flows. 



the DNS (see figure [5f| left), which indicates ODT to be an accurate near-wall 
model, yet by including additional 3D resolution iVLES within ODTLES the 
PDF is significantly improved. A more detailed investigation (see figure 
right) shows that ODTLES underestimates rare backflow events. The reason 
could be that the responsible 3D structures near the wall are not represented 
due to the coarse 3D resolution, because the 3D cell size in wall units is 
AxLES’+ 32 for the highest considered 3D resolution A^les = 64 (with 
Rcr = 1020). The result with lower Reynolds number [Rcr = 395) with 
3D cell size in wall units AxEes.+ ~ 24.7 supports this hypothesis. Here we 
as sume a low Reynolds num ber sensitivity of the wall shear stress statistics. 


as 


Schlatter and Orlii 


(2010) report. 


In summary ODTLES is able to capture the mean flow and turbulence 
statistics of highly turbulent flows up to Rcr = 10000 within a simple domain. 


The computational costs are significant lower compared to DNS: Lee et ah 


(2013) report to use about 260 x 10® CPU hours on 786 x 10^ cores for the 


DNS with Re-r = 5200, while the corresponding ODTLES simulation takes 
~ 10000 CPU hours on 24 cores. 

A detailed investigation of the expected computational costs relative to 
RaNS, LES and DNS in section shows ODTLES to be convenient to de¬ 
scribe highly turbulent flows in domains of moderate complexity. 

6.3. ODTLES: Square Duct Flow 

Sections 16. II and 16.21 show turbulent channel flow results for the ODTLES 
model to be in good agreement with DNS. With increasing 3D resolution 
a convergence to DNS is observed, but even within an ‘ODT-limit’ (corre¬ 
sponding to a single 3D cell) the ODT model reproduces key flow features 
stand-alone 


see 


Appendix A). 


In this section the square duct flow is investigated. This flow combines a 


simple geometry (see figure 6a) and a complex flow behavior including sec¬ 


ondary instabilities (secondary flow of Prandtl’s second kind) and turbulent 
fluctuations. These secondary instabilities represent a 3D flow phenomena 
which is not captured by the ODT model. For low Reynolds numbers (near 
the value for sustained turbulence) the size of the secondary flow structures 
in cross-stream direction corresponds to the half duct height (h). This flow 


25 



























regime in a square duct is investigated by Uhlmann et al. (2007) using DNS. 
Even for higher Reynolds number the cross-stream extension of the secondary 
instabilities is rather large scale compared to turbulent fluctuations occurring 
e.g. near the wall. Nevertheless these small scale fluctuations play an impor¬ 
tant role for the duct flow because they generate secondary instabilities. The 
ODT model was shown to accurately describe small scale fluctuations and 
the secondary instabilities correspond to a rather large scale 3D flow feature 
which can be described by the XLES framework. Thus ODTLES is a highly 


promising model to describe the duct flow behavior, as Gonzalez-Juez et al. 


(2011) and Glawe et al. (2013) showed in previous ODTLES studies. In these 
works no coupling terms (last line in Eq. (26)) between the XLES-grids are 
considered. 

The GEL number is chosen following Eq. (23) with CFL < 1. The 
ODT model parameters are C = 6.5, Z = 330, and which is 


identical to the channel flow setup in section 6.2 

The flow is averaged for taveUs/h > 2800 non-dimensional time units 
(with the bulk velocity ub) after reaching a steady state which is assumed to 
be sufficient to investigate the secondary instabilities. This is supported by a 
study of vortex structures by Uhlmann et al. (2007) where nearly symmetric 
8-vortex structures are observed for lower averaging times at the investigated 
Reynolds numbers Rcb > 2600. Additionally the ODTLES results are aver¬ 
aged over the 4 quadrants. Note that the DNS by |Pinelli et al. (2010) uses 
taveUB/h > 7000 to produce meaningful high order statistics. 

Figure [^compares the secondary flow computed by ODTLES for a moder¬ 
ate bulk Reynolds number Rcb ~ 2600 with the DNS by Pinelli et al. (2010) 
(online available: Uhlmann (2013)). Hereby ODTLES uses Ales = {16,32} 
3D large scale cells per direction and the XLES specific small scale properties 
are resolved by Arss = 512 cells. 

The primary flow and key features of the secondary flow are in good 
agreement with the DNS results even with the very low 3D resolution of 
Ales = 16 cells. Furthermore the ODTLES results indicate a convergence 
towards the DNS results with increasing 3D resolution. 

The secondary mean velocity field alters with increasing Reynolds num¬ 
ber, as investigated in figure In particular each secondary vortex structure 
tends towards a triangular shape for increasing Reynolds number. The vortic- 
ity approaches the duct corner with increasing Reynolds number, indicating 
fast changes in secondary flow directions in this area. 

ODTLES results show a tertiary instability above some threshold Reynolds 
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(a) Geometry and illustrated XLES- 
grids. 




(c) ODTLES: IVles = 16, Res = 2667. (d) ODTLES: IVles = 32, Res = 2514. 


Figure 6 : All results are averaged in time and streamwise direction. One 


quadrant of the duct is shown (see 6 a). The 3D grid is indicated by white 


lines (for the DNS only in the corner region). ODTLES properties are illus¬ 
trated like cell centered and show additional small scale features (resolved 
by AIrss = 512 cells) using the XLES-grid highly resolved in vertically X 2 - 
direction (horizontally xa-direction) in the lower right (upper left) triangular 
region, as illustrated in 6 a Contour lines of the primary mean flow (black) for 
Ml = {0.2, 0.4, 0.6} max(Mi), streamlines of the secondary mean flow (^ 2 ,^ 3 ) 


in red and the 2D vorticity U 2 d 
shown. 


= — dx^U 2 (RGB color coded) are 
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(a) DNS: 

Res = 3500. 



(c) ODTLES: 

Rgb = 8446, iVj^gs = 2048. 



(b) ODTLES: 

Ree = 3485, iVRss = 1024. 



(d) ODTLES: 

Res = 17338, Nrss = 2048. 


Figure 7: ODTLES and DNS square duct results with iVLES = 32 for different 
Reynolds numbers Res- Primary and secondary streamlines and vorticity are 
illustrated similar to figure 






















































(a) Res = 3485, N^ss = 1024. (b) Res = 17338, JVrss = 2048. 


Figure 8: Zoom into the corner region of figure 7b {Res = 3485) and 7d 
{Rcb = 17338). Primary and secondary streamlines and vorticity are il¬ 
lustrated similar to figure ODTLES properties are illustrated like cell 
centered (this leads to the gap in the flow held at the corner). 8a addi¬ 
tionally shows the small scale resolution (illustrative with A^rss = 256). An 
exemplary cross-stream how approaches the corner parallel to the horizon¬ 
tal wall (green arrows) with wall-normal velocity gradients highly resolved 
only in the XLES-grid 2 (compare to hgure [6a]) . Near the corner the vertical 
wall forces a how stagnation and an associated pressure gradient drives the 
how in the horizontal direction (blue arrows) which is highly resolved only 
in XLES-grid 3. For high Reynolds number duct hows (8b) both ehects oc¬ 
cur within one 3D cell, which is not well represented by ODTLES because 
the coupling procedure only communicates small scale ehects ahecting the 
large scale. This could lead to unphysical how behavior within the 3D cell 
containing the corner. 
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number (between 8446 and 17338) in the corner region. This model result is 
not conclusive because the tertiary instability is not resolved sufficiently in 
3D (see £gure[^. 

The 3D tertiary structure in the corner region is of similar size as the 
corresponding 3D computational cell. Its influence on the flow is primarily 
local. The insufficient 3D resolution possibly leads to an under-resolved 
pressure gradient which could prevent an adequate change of a fluid parcel’s 
direction (blue arrows within £gure|^. Nevertheless the cross-stream flow is 
adequately resolved in the wall-normal direction (green arrows within hgure 
1^. More reliable and conclusive results might be possible through future 
DNS studies. 

ODTLES demonstrated its ability to describe and predict non-trivial flow 
behavior including secondary instabilities within a duct flow. We suspect that 
the flow transition leading to tertiary instability is a model artifact, although 
a physical cause cannot be ruled out until dehnitive evidence such as a DNS 
result becomes available. Nevertheless the tertiary structure is very local and 
not preventing the ODTLES model from describing the key flow features of 
the primary and secondary flow. 


7. ODTLES: Resolution Properties and Efficiency 


Different turbulence models, e.g. RaNS, wall-modeled LES (LES includ¬ 
ing a near-wall model), wall-resolved LES (LES with near wall resolution) 
and ODTLES, differ strongly in both represented physical effects and compu¬ 
tational effort. In this section the computational costs of the different model 
approaches are estimated by developing a relation between the grid-size (used 
as a measure for the computational effort) and the Reynolds number follow¬ 


ing Chapman (1979) and especially Choi and Moin (2012) and references 


cited therein. 

The investigated domain is a box of size Li x L2 x L^. A highly turbulent 
boundary layer over a flat-plate airfoil of the thickness 6 fills the volume 
[xq, Li] X 5(a;i) x L 3 . The flow is assumed to reach the plate at Xi = Xq with 
a turbulent intensity The boundary layer size 5(xi) increases until 

reaching xi = Li with a corresponding Reynolds number Re^^. The number 
of grid cells N within the turbulent boundary layer is estimated for ODTLES 
and compared to RaNS, LES and DNS. 


From Choi and Moin (2012) we extract the Reynolds dependent grid size 
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for RaNS and wall-modeled LES: 


-^3 2/7 

NRaNS\wm = —nin2n^Re J 
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(43) 


for a wall-resolved LES 


AT n noi -^3 TT'2,laminar j-. 13/7 

iV^. = 0.021- 
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(44) 


and for DNS 


Ndns = 0.000153^i?ef^^ 

-Li 


Re 


xo 


Rei, 


23/14 


(45) 


Here n^nyTiz is the nnmber of grid points within the cnbe 6 x S x 6, n 2 ^iaminar 
is the number of wall-normal grid points within the laminar sublayer, and 


is the LES cell size in wall units. 


Following Chapman (1979) RaNS typically resolves the cube 6x6x6 
using n^nyUz ~ 1 x 20 x 0.5 = 10 cells, while for wall-modeled LES [Choi and 


Moin (2012) report typical grid resolutions n^UyUz ~ [1200,33000]. In wall- 


resolved LES 


Choi and Moin 


(2012) find typical resolution values 


n2,i 


aminar 


AxJ 




[— -1 

Here LES models represent turbulent scales down to the inertial range of 
the turbulent cascade. The ODT model, applied within XLES, potentially 
describes the full turbulent cascade within a ID sub-domain, which leaves 
the 3D grid to capture non-turbulent effects (e.g. the domain or secondary 
instabilities). For the flat-plate airfoil even ODT stand-alone potentially 


leads to reasonable results for the case of a turbulent boundary layer (Lignell 


et ah (2013) apply ODT to a comparable turbulent case, but including buoy¬ 


ancy, by spatially advancing the ODT line). In consequence the (equidis¬ 
tant) XLES 3D resolution N^. in x^-direction is chosen independently of the 
Reynolds number (unless Reynolds-number variations triggers a global flow 
structure transition, like the secondary instabilities in a turbulent duct). 

In the current ODTLES implementation the resolved small scales are 
represented by iV^sSk equidistant cells in Xfc-direction {k = {1,2,3}). Fol¬ 


lowing Choi and Moin (2012) for highly turbulent flows the number of grid 
points resolving the Kolmogorov length scale along a small distance dxi is 
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Ni = Equidistant ODTLES uses the smallest length scale 

X\ 

globally in all ID sub-domains, leading to 


Nqdtles — 


(46) 


A factor Kqdt ~ 3 x 6 takes into account that 3 XLES-grids are used (we 
assume that Ni = N 2 = A 3 ) and equidistant ODT uses at least 6 cells to 
allow a turbulent event (eddy) within the Kolmogorov scale. 

In principle ODTLES can be extended to non-equidistant grids within the 
ID sub-domain, which is for example realized by the adaptive ODT (aODT) 
implementation by Lignell et ah (2013). Although adaptive ODT is not 
used as a sub-grid model within an XLES approach yet, we investigate this 
interesting case as a worthwhile perspective and refer to it as aODTLES. For 
an adaptive grid we assume on average a resolution similar to DNS (in ID) 
and integrate over the boundary layer thickness with - = 0.16i?ea; ' (see 


Choi and Moin (2012)) in the ID sub-domain, leading to 

Re 


NaODTLES — 0.0103936AaODT Ai A 3 — 


1 - 


XO 


Re 


Li 


25/14 


(47) 


For adaptive ODT, note that we assume KaODT = 3 because 3 XLES-grids 
are required (here we assume each XLES-grid uses the same RSS resolution 
and Ai = A2 = A3). There is no additional assumption of a minimum 
number of cells representing the Kolmogorov length. 


We compare typical RaNS and LES resolutions (following Choi and Moin 


(2012) and Chapman ( 1979[ )) with the ODTLES and aODTLES approach 
for different 3D resolutions in hgure Additionally numerical computations 
for a turbulent channel are shown assuming a similarity of the turbulent 
Reynolds number in the channel and Rei^- Hereby the DNS by Lee and 
Moser (2014b) (we assume A ^ 8.5 x 10®, estimated for L 3 /L 1 = 4), an 


aODT result by Meiselbach (2015b) {NaODT ~ [80000,120000] (Meiselbach 


(2015a))) and the ODTLES result with Re^ = 10000 in section 6.2 are used. 

For weakly turbulent flows, ODTLES is subject to additional computa¬ 
tional costs compared to standard LES. But ODTLES requires 3D resolution 
independent of the turbulence intensity (except secondary effects), and thus 
highly turbulent flows in moderately complex domains are well described 
with low computational costs. In some flow regimes ODTLES is more effi¬ 
cient than wall-modeled LES and although it represents advective and dif¬ 
fusive effects down to the Kolmogorov length scale. Incorporating adaptive 
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Figure 9: Number of grid points N required for numerical simulations of a 
flow over a flat-plate airfoil with aspect ratio Fs/Li = 4 and a turbulent in¬ 
flow with Rcxq = 5 X 10®. Wall-modeled LES (^ 1712^3 = 2500) , wall-resolved 
LES {n2,iaminar/ = 1/200), RaNS (nin 2 U 3 = 10), ODTLES and 
aODTLES (ODTLES with adaptive ODT) with Ni = = {10,100} cells, 

and the ODT and adaptive ODT (aODT) stand-alone model. Additionally, 
actual simulation cases for a turbulent channel (assuming Returb ~ Rli) are 
shown for DNS (Rcr = 5200, black point), ODTLES {Rcr = 10000, green 
point), and aODT {Rcr = 6 x 10®, magenta point). 
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ODT into a XLES framework seems to be an especially promising alternative 
to wall-modeled LES and even RaNS simulations for highly turbulent flows 
in simple domains (requiring low 3D resolutions). 

Note that in the presented estimation the costs of the computation within 
one discrete cell is neglected, because the different modeling strategies typi¬ 
cally vary by a low factor (< 2), which is not strongly affecting the estimation 
in figure 


8. Conclusions 


Part I (Glawe et ah (2015)) introduces XLES, an extended LES approach, 
which is a new strategy to simulate complex turbulent flows. The ODTLES 
model is one special approach in the XLES family of models, employing ODT 
as a sub-grid model. XLES in general and especially ODTLES are designed 
to describe highly turbulent flows in domains of moderate complexity. These 
problems especially occur in fundamental research studies of e.g. atmospheric 
flows and are also relevant in engineering. 


A previous ODTLES version (by Schmidt et ah (2008) and Gonzalez-Juez 


et ah (2011)) can also be interpreted as ODT closed XLES, but with XLES 
velocities only coupled by a pressure projection. Thus oscillations in root 


mean square velocities occur (reported by Gonzalez-Juez et ah (2011)), but 
without significantly affecting mean profiles. The introduced XLES coupling 
terms directly carry over to scalar properties residing on different XLES-grids 
(e.g. required for a heated duct flow) which is a unsolved problem within the 
previous ODTLES formulation. 

XLES time advances multiple coupled 2D hltered Navier-Stokes realiza¬ 
tions, each having one Gartesian direction that is highly resolved. A one¬ 
dimensional modeling approach like ODT takes advantage of the specihc 
symmetry of a 2D filter. 

Especially within ODTLES, the ability of ODT to describe the full tur¬ 
bulent spectrum allows strongly reduced 3D resolutions without corrupting 
key flow features. 

In XLES, separated physical effects (contrary to separated 3D scales in 
LES) are represented by appropriate approaches: the 3D resolution repre¬ 
sents the 3D domain and other physical effects not captured by ODT, e.g. 


secondary instabilities (see section 6.3), while turbulent effects, not captured 


by XLES, are appropriately represented by the ODT sub-grid model. This 
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includes the representation of molecular diffusion and turbulent advection at 
the Kolmogorov length scale (within a ID sub-domain). 

ODTLES accurately describes a turbulent channel flow up to friction 
Reynolds number Re < 10000 with high accuracy, even with coarse 3D res¬ 
olution (e.g. IVles = 16 cells per direction) and is able to reproduce the 
primary and secondary flow in a square duct with similar 3D resolution. 

The focus of this work is to introduce the mathematical framework neces¬ 
sary to derive the XLES approach, which is one possible way to incorporate 
ODT into 3D simulations. This is implemented here for a simple ODT model 
to avoid e.g. interpolation effects that are introduced by the adaptive ODT 


model by Lignell et ah (2013). Adaptive ODT outperforms the turbulence in¬ 


tensity reachable by equidistant ODT, as Meiselbach (2015b) recently showed 


using adaptive ODT simulations up to Rcr < 6 x 10^. 

Incorporating this adaptive ODT model into the XLES hlter approach 
can potentially serve as an alternative to RaNS simulations in industrial 
applications and additionally include a wide range of small scale physical 
effects, e.g. for flows including Prandtl number effects and combustion. 

Various additional physical effects, e.g. additional scalar helds, are well 
tested for ODT and can easily be adapted to ODTLES. Furthermore the 
XLES framework is the hrst approach which consistently couples scalar prop¬ 
erties and velocities between XLES-grids (in the previous ODTLES version 
only velocities were coupled). This enables ODTLES to compute e.g. fun¬ 
damental meteorological flows. 

Thus ODTLES accurately describes highly turbulent flows including fully 
resolved small scale effects with low computational costs for simple domains. 
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Appendix A. ODT: Illustrative Results 


As discussed in section]^ ODT error terms ctodt can roughly be estimated 
by comparing ODT and DNS flow statistics. The turbulent channel is an 
appropriate study case for ODT because of its distinct predominant direction. 
ODT results with friction Reynolds number Rcr = 395 are compared to 
DNS (Kawamura (2014)), which allows additional comparison to LES-U and 
XLES-U results in (Glawe et ah, 2015, section 3.4) and ODTLES results in 
section 16.1[ 


The ODT model parameters are C = 6.5 and Z = 300 which are the 
same as the ODTLES parameters in section 6T and 6^2 The maximum 
eddy length is chosen to equal the channel half height h. The ODT 
resolution is Nqdt = 1024. To compute reliable ODT flow statistics a larger 
averaging period (or ensemble averaging) is required: the average time is 
tave = 12800 non-dimensional time units after reaching a steady state, which 
is significantly larger than tave = 20 in DNS. 

Figure A.10| illustrates some representative results for ODT; The averaged 
streamwise velocity profile (see figure A. 10a) and the overall turbulent kinetic 
energy (see figure A.lOc-A.lOe) are described very well within the ODT 
model including the full spectrum of the turbulent cascade. The similarity of 
ODT and DNS results implies a considerably small ODT model error uodt- 
The ODT computing time is only 9 CPU-seconds to simulate for tave = 20 
(the overall computing time is higher, because in the example the flow is 
averaged over tave = 12800). ODT is a convenient sub-grid model because 
of its low computational costs. Thereby ODT can compute a wide range of 
complex physical effects. In the turbulent channel case, ODT dynamically 
produces realistic wall profiles (see section 4.2). These properties outperform 
commonly used eddy viscosity models. 

Further ODT results, including various physical small scale effects, are 


available in the literature (e.g. see 

Kerstein et al. 


>001), 

Schmidt et al. 

(2013) 

Wunsch and Kerstein (2005) 

, and 

Schulz et al. 

(2013 

)• 


Appendix B. ODTLES: Numerical Implementation 


Table B.l| summarizes the numerical discretizations of the single frac¬ 
tional steps within the ODTLES advancement cycle (section]^. The ODT 
advancement (Eq. (p^) contains the XLES diffusion terms (^xles) treated 
as forcing terms for ODTLES. 
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(a) ODT: (b) Streamwise (wrms) and spanwise 

law of the wall. ('U^rms) velocity RMSs. 




(c) Production (Prod) and Dissipation (d) Viscous transport of the turbulent 
(Diss) of the turbulent kinetic energy, kinetic energy (tv). 



(e) Advective transport of the turbulent 
kinetic energy (ta). 

Figure A. 10: Turbulent channel flow results [Rct = 395) for DNS (small 
crosses) and ODT (solid) with Nqdt = 1024 cells. 
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Table B.l: Numerical Schemes: 

In time: EEl (1st order explicit Euler), RK3 (3rd order Runge-Kutta), CN 
(2nd order Crank-Nicolson), lEl (1st order implicit Euler). 

In space: UPl (1st order upwind), CDM (2nd order central difference 
method). 


factional step 

time scheme 

spatial scheme 

Eq. ( 

36 

) 

EEl 

UPl 

Eq. ( 

37 

) 

RK3-CN + RK3-RK3 

CDM 

Eq. ( 

38 

) 

EEl 

CDM 

Eq. ( 

39 

) 

lEl -|- triplet map -|- EEl 

CDM 

Eq. ( 

40 

) 

EEl 

CDM 


The coupled advection schemes RK3-CN-CDM and RK3-RK3-CDM are 
described and validated in part I (Glawe et ah, 2015, section 3.3). All veloc¬ 


ities are discretized using a staggered grid. 
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